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Abstract 



1—5 I We study microwave background anisotropies induced by scaling seed pertur- 

00 ' bations in a universe dominated by cold dark matter. Using a gauge invariant 

^N , linear perturbation analysis, we solve the perturbation equations on super- 
horizon scales, for CMB anisotropies triggered by generic gravitational seeds. 

>- ■ We find that perturbations induced by seeds — under very mild restrictions — 

^^ . are nearly isocurvature. Thus, compensation, which is mainly the consequence 

f"^ I of physically sensible initial conditions, is very generic. 

CN ' We then restrict our study to the case of scaling sources, motivated by 

^^ ' global scalar fields. We parameterize the energy momentum tensor of the 

Q"^ . source by "seed functions" and calculate the Sachs- Wolfe and acoustic contri- 

fH I butions to the CMB anisotropies. We discuss the dependence of the anisotropy 

Q^' spectrum on the parameters of the model considered. Even within the re- 

Q . stricted class of models investigated in this work, we find a surprising variety 

|H I of results for the position and height of the first acoustic peak as well as for 

52 I the overall amplitude. In particular, for certain choices of parameters, the 

IIV ' spectrum resembles very much the well known adiabatic infiationary spec- 

• i-H . trum, whereas for others, the position of the first acoustic peak is significantly 

rS ' shifted towards smaller angular scales. 



PACS numbers: 98.80-k 98.80.Hw 98.80C 



1 Introduction 

The origin of the large scale structure in the universe, is clearly one of the most 
important open questions in cosmology. Within the framework of gravitational in- 
stability, there are two currently investigated families of models to explain the for- 
mation of the observed structure. Initial density perturbations can either be due to 
"freezing in" of quantum fluctuations of a scalar field during an inflationary period 
|I|], or they may be seeded by topological defects, which can form naturally during 



a symmetry breaking phase transition in the early universe 0. Inflationary fluc- 
tuations are produced at a very early stage of the evolution of the universe, and 
are driven far beyond the Hubble radius by inflationary expansion. Subsequently, 
they are not altered anymore and evolve freely according to homogeneous linear 
perturbation equations until late times. These fluctuations are termed "passive" 
and "coherent" |Q. "Passive", since no new perturbations are created after infla- 
tion; "coherent" since randomness only enters the creation of perturbations during 
inflation, subsequently they evolve in a deterministic and coherent manner. 

On the other hand, in models with topological defects or other types of seeds, 
fluctuations are generated continuously and evolve according to inhomogeneous lin- 
ear perturbation equations. The seeds are any non-uniformly distributed form of 
energy, which contributes only a small fraction to the total energy density of the 
universe and which interacts with the cosmic fluid only gravitationally. We will be 
particularly interested in the case of global topological defects, playing the role of 
seeds. The energy momentum tensor of the seed is determined by the defect (seed) 
evolution which, in general, is a non-linear process. These perturbations are called 
"active" and "incoherent" 0]. "Active" since new fluid perturbations are induced 
continuously due to the presence of the seeds; "incoherent" since the randomness 
of the non-linear seed evolution which sources the perturbations can destroy the 
coherence of fluctuations in the cosmic fluid. 

The cosmic microwave background (CMB) anisotropics provide a link between 
theoretical predictions and observational data, which may allow us to distinguish 
between inflationary models and defect scenarios, by purely linear analysis. On 
large angular scales, both families of models predict an approximately scale-invariant 
Harrison-Zel'dovich spectrum [|, Q. Although, perturbations from defect models 
are non-Gaussian, this signature is probably rather weak, especially on large scales, 
where cosmic variance is substantial, and its observation might be quite difficult. 

Acoustic peaks have been extensively studied in inflationary models, where obser- 
vations of amplitude and position of the peaks can be used to determine cosmological 
parameters [Q. Some studies of simplified models where perturbations are seeded 
by topological defects, have already appeared in the literature]^, |^. 

In this paper, we present a general investigation of acoustic peaks for models 
with active perturbations. We estimate the Sachs- Wolfe and acoustic contributions. 
Our "seed functions" , which determine the energy momentum tensor of the source 
are motivated from 3d numerical simulations of ti^ defects, textures |P, 0, in a 
universe dominated by cold dark matter (CDM). We restrict ourselves to scalar 
perturbations. 

In section 2, we study CMB anisotropics triggered by generic gravitational seeds. 
In particular, we present the equation for the coefficients Ci and discuss the Sachs- 
Wolfe contribution, the acoustic peaks and Silk damping. We solve the perturbation 
equations for super- horizon scales. Studying the Bardeen potentials \E' and $, which 
describe the scalar geometry perturbations, we find that compensation is automat- 
ically obtained, i.e., that perturbations are nearly isocurvature. In section 3, we 
restrict our investigation to the case of global scalar fields, for which we deduce 
the power spectra of the seed functions from numerical simulations and analytical 
scaling arguments. The results obtained there apply for general scaling sources. In 
section 4, we present some numerical examples and discuss how the characteristics of 
the acoustic peaks depend on the model. We summarize our conclusions in section 
5. . 

Notation: The Friedmann metric is given by a?{—dt^ + ^jijdx^dx^), where a 

denotes the scale factor, t is conformal time and 7 is the metric of a three space 
with constant curvature K. We shall consider a universe dominated by cold dark 
matter and discuss the case K = exclusively. An over-dot stands for derivative 
with respect to conformal time t, while prime denotes the derivative with respect to 
kt = X. Subscripts (or superscripts) ,, {^^^) and c (^'^■') indicate the baryon-radiation 



plasma and CDM, respectively. 



2 CMB anisotropies from seeds 



The coefficients C^ represent the angular power spectrum of CMB anisotropies. They 
can be given in terms of the expansion of the angular correlation function 

(f(n)f(n')) I <„..,_., = i^E(2^+ l)^^^Kcos^) . 

We want to investigate these coefficients in models where the fluctuations are induced 
by seeds. We restrict ourselves to scalar perturbations, but this analysis is easily 
extended to include vector and tensor contributions. 

If we neglect Silk damping in a ffist step and integrate the photon geodesies in 
the perturbed metric, gauge invariant linear perturbation analysis leads to [ITT], IT2 
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-iDM(x)-l/.(xy-(v[;-$)(x)'^ 
+ (\i>-^){^,t')dr, (1) 



where over-dot denotes derivative with respect to conformal time t. $ and \l/ are the 
Bardeen potentials, quantities describing the perturbations in the geometry, V is the 
peculiar velocity of the baryon fluid with respect to the overall Friedmann expansion 
and D^^^ specifies the intrinsic density fluctuation in the radiation fluid. There are 

several gauge invariant variables which describe density fluctuations; they all differ 
substantially on super-horizon scales but coincide inside the horizon. Dg corresponds 
to the density fluctuation in the so-called "flat slicing" , where the perturbation of the 
3-dimensional Riemann scalar vanishes. The initial time is the time of decoupling, 
tdec of baryons and radiation, which occured at a redshift of Zdec ~ 1100. 

The final time values in the square bracket of Eq. (|ip give rise only to monopole 
contributions and the dipole due to our motion with respect to the CMB and are 
disregarded in the following. Taking the Fourier transform of Eq. (|l|), we then 
obtain 

^(k,n) = e^('^")*"[iDM(k, W) + V(k, W) • n+ (^ - $)(k, W) 



T ^ ' ^ 'A 



+ r°(^-<i>)(k,t)e-'(^-")*cit] . (2) 



The ffist term in Eq. @ describes the intrinsic inhomogeneities on the surface of 
the last scattering due to acoustic oscillations prior to decoupling. It also contains 
contributions to the geometrical perturbations 0. The second term describes the 
relative motions of emitter and observer. This is the Doppler contribution to the 
CMB anisotropies. It appears on the same angular scale as the acoustic term and we 
denote the sum of the acoustic and Doppler contributions by "acoustic peaks" . The 
last two terms are due to the inhomogeneities in the spacetime geometry; the ffist 
contribution determines the change in the photon energy due to the difference of the 
gravitational potential at the position of emitter and observer. Together with the 

part contained in D^^'' they represent the "ordinary" Sachs- Wolfe effect. The second 



term accounts for red-shifting or blue-shifting caused by the time dependence of 
the gravitational field along the path of the photon (Integrated Sachs- Wolfe (ISW) 
effect). The sum of the two terms is the full Sachs- Wolfe contribution (SW). 

On angular scales 0.1° ~ ^ ~ 2°, the main contribution to the CMB anisotropics 
comes from the acoustic peaks, while the SW effect is dominant on large angular 
scales. For topological defects, the gravitational contribution is mainly due to the 
ISW. The "ordinary" Sachs Wolfe term even has the wrong spectrum, a white noise 
spectrum instead of Harrison-Zel'dovich [^ . 



From Eq. (||) the C^'s are found to be 
<|A,(k)P>,_, 



^'=-JikW'- ''' 



with 



2i + l 



Je{kto) 



^L>f (k, W) + (^ - $)(k, W)] - j;(A;to)V,,(k, W) 
+ r {i> ~^){\<i,t')n{k{t,-t'))dt' 

''i'dac 

-Dj^)(k, W)j£(Ho) - j;(A;to)V,(k, Uec) 

+k r (^ - ^)ik,t')fMto-t'))dt' ; (4) 



ji denotes the spherical Bessel function of order i and j'^ stands for its derivative 
with respect to the argument. 

On scales smaller than about 0.1°, the anisotropics are damped due to the fi- 
nite thickness of the recombination shell, as well as by photon diffusion during re- 
combination (Silk damping). Baryons and photons are very tightly coupled before 
recombination and oscillate as one component fluid. During the process of decou- 
pling, photons slowly diffuse out of over-dense into under-dense regions. To fully 
account for this process, one has to solve the Boltzmann equation (see, e.g. []T2|). A 



reasonable approximation can however be achieved by multiplying the A^ with an 
exponential damping envelope 'D^{k) which is given in Ref. ||13| . 

We now discuss the calculation of the C^'s for perturbations with seeds. Since 
the background contribution of the energy momentum tensor of the seeds vanishes, 
its components Q^u are gauge invariant perturbation variables. They can be de- 
composed into scalar, vector and tensor contributions. Here we restrict ourselves 
to scalar perturbations. We express the scalar degrees of freedom of Q^v in terms 
of the gauge invariant perturbation variables fp,fp,fv,fn, which parameterize the 
energy density, pressure, scalar velocity potential and anisotropic stress potential of 



seeds, respectively (see [|TT], |T2| 

Ooo = MVp (5) 



©1? = M'U, (6) 

eif = M2[{/,-(l/3)A/j7., + /.,d. (7) 

where A denotes the Laplacian with respect to the metric 7 of the three space and 
M is a typical "mass", energy scale, of the seeds. The superscript ^^^ indicates 
that only the scalar contribution to Ojo and 0,^ is obtained in this way. Numerical 
simulations show that the vector and tensor perturbations make up about 20% of 



the energy moinentuin tensor on super-horizon scales |Ty]. Since we assume that the 
seeds interact with other matter components only gravitationally, the seed functions 
satisfy the following covariant conservation equations [|12| 



/,-A/„ + (a/a)(/, + 3/p) = 
/, + 2(d/a)/,-/p-(2/3)AA = 0. 



(8) 
(9) 



We consider the matter content of the universe as a two-fluid system: the 
baryons+radiation plasma, which prior to recombination is tightly coupled, and 
cold dark matter (CDM). Before recombination, the evolution of the perturbation 
variables in a spatially fiat background, f2 = 1, is described by [|14 
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a ^ 

Df^ + kV, 

^2 



K + (1 - 2,cl)-Vr - k{m - 3c^$) - k- 



1 + Wr 



-Pt^ 



Vc + -K - km 







0, 



(10) 

(11) 

(12) 
(13) 



where subscripts r-, c (superscript ^''\ '^^^) denote the baryon-radiation plasma and 
CDM, respectively; D and V are density and velocity perturbations; w = Pr/Pr, 
c1 = Pr/pr and p = Pr + Pc- The geometrical perturbations \E' and $ can be separated 



into a part coming from standard matter and radiation (subscript 
due to the seeds (subscript s)- 



and a part 



^ = *„ + ^. (14) 

$ = $„ + $., (15) 

where \E's and <l>s are determined by the energy momentum tensor of the seeds. We 
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$. 



^r 



^^[PrOp + p^Df - 3{pril + m,) + pj* 
a 



+3-k-^{pr{l + Wr)Vr + PcVc}] 



-$. 






-$, - 2e/^ , 



(16) 
(17) 
(18) 
(19) 



where e = AirGJVP. This parameter has to be small to validate cosmological pertur- 
bation theory. In other words, the mass M has to be significantly smaller than the 
Planck mass. For global scalar fields, it actually turns out that the typical amplitude 
of geometrical perturbations is of the order of e, so that the COBE normalization 
requires M ~ lO^^GeV. In this work, we neglect the contribution of neutrino fluc- 
tuations. Anisotropic stresses in the matter components are explicitly set to zero, 
Ilm = 0, which implies Eq. (p!7|). The anisotropic stresses in the source, f.„, can 
therefore not be compensated. 



To solve the above system of equations, we need to specify initial conditions. For 
a given scale k, we choose the initial time tin early enough, such that the perturba- 
tions are super-horizon and the universe is radiation dominated at tj„. We set x = kt 
and denote by a prime the derivative w.r.t x. The super-horizon limit is thus the 
hmit X <^ 1. Choosing {x,k) as independent variables, the perturbation equations 
reduce in this limit to 

Df + ^K = (20) 

Df + V, = (21) 

V; + (^ - ^) ~ ^Dl[^ = (22) 

V^ + Vjx-^ = (23) 

^m + ^m = (24) 

$. = ek-\f, + 3fjt) (25) 

2 

$ = ^$, + V;./x (26) 

^ + $ = 2eU. (27) 

If fp, /„ and f-jr are different iable in the vicinity of a; = 0, we can solve the above 
system exactly. We take the derivative of Eq. ( ^0]) and replace V^' with Eq. (^2]) . 



Using Eqs. (|6|, |7|, ||) we then find 



oi;y' + lD';y = .(2f, + ^f, + l^). (28) 

Differentiability now guarantees that the source term on the right hand side of Eq. 
( pSD is given by eA{k)x°' + higher orders. We then obtain to lowest order in x 

dM = ^^(^^ a;"+' f29) 

_ 3eA{k) -^ 
^^ " A{a + 3f ^^^> 

^-^ - ^!4M^.«. (31) 

4(a + 3) ^ ' 

On the other hand, the seed perturbations are of the order of 

$, , ^, oc eA(A;)x°-2 > $ , ^ , (32) 

if /,r ~ {x^ /'ik'^)fp + {x/k)fy. For scaling sources, we shall see that these two terms 
are of the same order of magnitude. In other words, $ <^ $s and, if /^r is not 
extremely large, \1/ -C \E's on super-horizon scales. The main reason for this find- 
ing certainly lies in choosing the correct initial conditions which have to vanish in 
the absence of sources. We could always add a homogeneous contribution to D^^^ 
which would destroy this behaviour. We consider this choice of initial conditions 
as the most natural way to obtain compensation: the presence of matter and radi- 
ation reduces the Bardeen potentials on super- horizon scales by a factor x"^. Only 
the contribution $ -|- \l/ = — 2e/,r, which is due to anisotropic stresses, cannot be 
compensated by matter and radiation. If /^ 7^ 0, there is compensation provided 
/^ < max(/pt2 , fj). 



We therefore conclude that seeds, which are uncorrelated on super-horizon scales, 
are compensated by the presence of matter and radiation, where we define compen- 
sation as the suppression of the total Bardeen potentials by a factor x^ with respect 
to <l's,^s- In this sense, the type of seed perturbations discussed here are nearly 
isocurvature fluctuations. 

Within the context of scaling sources (seed functions with white noise spectra), 

the basic ingredient which leads to compensation, is not the absence of perturbations 
on very large scales or on very early times, but the fact that we only consider the 
particular solution of the second order differential equation for the perturbation 
variable D'^''^ . Clearly, a homogeneous contribution to D'^^^ can destroy this finding. 
In our case, perturbations are induced by the presence of the seeds, and therefore 
these initial conditions are the most physical ones. 

This result is important, since compensation has usually been understood either 
as a consequence of the integral constraint |l^, [1^, |^ or as a consequence of causality 
of the source perturbations [p!7|] . In our work compensation arises naturally for scaling 



sources and it can be generalized to sources with arbitrary spectra which satisfy 

U < maxifpt' , U) . (33) 

Clearly for 11 = 0, which is the most natural assumption for non-relativistic cosmic 
fluids, /tt has to be small since it cannot be compensated, $, \E' ~ (9(e/^). If f,, is 
larger than the limit given in Eq. (BSp, then 11 7^ is a necessary but not sufficient 



condition for compensation to occur. We know of one example, namely relativistic 
collisionless particles, where compensation can take place for certain choices of /^, 
due to the presence of anisotropic stresses 11 (see appendix). Collisionless particles 
are special in that they interact with each other and with the seeds only through 
gravity. In general, if particle interactions other than gravity determine 11, we do 
not expect compensation, since by definition the seeds interact with the cosmic fiuid 
only through gravity. 



3 CMB anisotropies induced by scaling sources 



We now restrict our study to scaling seeds. We first discuss as motivating example 
global scalar fields which, depending on their number of degrees of freedom, can 
lead to global topological defects during a symmetry breaking phase transition in 
the early universe . We shall, however, only make use of the general behaviour of 
the seed functions /, which we call "scaling" . In the absence of any intrinsic length 
scale other than the cosmic horizon, this is the behaviour which the seed functions 
assume by dimensional reasons. 

We consider an iV-component scalar field with potential V = A(0^ —jfY. In the 
(T-model approximation, the equation of motion for can be expressed solely in the 



dimensionless variable [3 = (p/i] ||T0| . In terms of [3 the energy momentum tensor of 
the scalar field is given by 



.2 In n 1 .,A 



Q,. = v\f^,,l3,.-7^9,Ax) ■ (34) 

Defining M^ = rf, the functions /, result in 

/p(k) = \:F[^' + {V(3f 



p'-\iyp?) 



/p(k) = \:f 



j-kv^) — 7m ii*^ 
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(3nP.J-\5^,iyPf 



(35) 



where J^[g\ denotes the Fourier transform oi g^ defined by g{k) = V^^l"^ J e^^^g{x)d^x. 
In what follows, the Fourier transform of the seed functions ^[/.j, will be denoted 
simply by /,. 

On super-horizon scales (3 and (3 are assigned random initial values, so they have 
white noise spectra initially. V/5 clearly has a fc^-spectrum. However, using the 
convolution theorem one finds that (V/9)^ has a white noise spectrum. Therefore, 
both fp and fp have white noise spectra on super-horizon scales. From the above 
expressions for /^ we find 



/. = ^^^/c^^g/3(q)/3(q-k)(g- A:), 



Expanding this expression in lowest order in k using that (3 and (3 have white 
noise spectra, we find that the term of order kP in the integral vanishes and the 
lowest order contribution to the integral is linear in kj, so that /^ also has a white 
noise spectrum on super-horizon scales. By similar arguments one can deduce that 
/^ has a white noise spectrum on super-horizon scales. 

The dimensions of fp and fp in physical space are (length)"^, therefore in k-space, 

fp and fp have dimensions (length)^^/^. Since on super-horizon scales {kt <^ 1), these 

functions have white noise spectra, they must behave as l/\/t. The corresponding 

arguments lead to a super-horizon behaviour for /^ oc ^/t and /tt oc t^'"^. So the 
power spectra of the seed functions behave like 



(/pp) = 


= Aj t-^ Fi{x 


{U?) - 


= Aj t-' F2{x 


(I/.') = 


= AltFsix) 


{U') - 


= Alt'F,{x) 



(36) 

where we choose the dimensionless constants Ai to be positive and such that -Fi(O) = 
1. The power spectra of the functions /, do not depend on the direction of k, thus 
the Fj's are even functions oi x = kt. Furthermore, since the energy momentum 
tensor of the source decays inside the horizon, we know that Fi ^ for x — ;> cx). 
This behaviour of fp, fp and /„ has also been found by numerical simulations. 

The temporal behaviour of fp and fp can also be understood from the following 
argument: the k = component of fp just corresponds to the average energy density 

multiplied by y/V and is thus proportional to V^^'^/t'^. On super-horizon scales, fp{k) 
is white noise superimposed on this average. The number N of independent patches 
in V is V/t^ and hence the amplitude of fp{k) is proportional to V~^^'^/{t'^N~^^'^) oc 
^-1/2 rjj^g same arguments hold for fp. From numerical simulations |T^ for tt^ 

defects, global textures, one finds that the average of /3^ over a shell of radius k can 
be modeled on super-horizon scales by 

(|/?T)(fc,t)~^. (37) 



We define a seed to be scaling if tlie power spectra of the seed functions behave 
as in Eq. (^61). We expect this scaling behaviour to be vahd not only for global scalar 
fields, but also for (local) cosmic strings. However, since the only decay mechanism 
for cosmic strings is through emission of gravitational radiation, we expect the func- 
tions Fi to decay slower on sub-horizon scales, than in the case of global fields, which 
decay very efficiently into Goldstone bosons. 

The system given by Eq. pO| ) to (pTf ) can be solved analytically if the stochastic 
variables /, are replaced by the square root of their power spectra. The results 
are thus to be taken with a grain of salt. But we believe that the r.h.s. of the 
following equations are good approximations to the square roots of the power spectra 
of the corresponding stochastic variables on the l.h.s., since, as we argue below, this 
coherence assumption does not significantly infiuence the results on super-horizon 



scales. Inserting the square roots of Eqs. (|36D in the system (^) to (pT]), one finds 
in the limit x <^1 

$, = e(Ai + 3A3)A;-='/'x-^/2 ^gg^ 

D^P = ^ei2A, + U, + A,)k~'/^xy' (40) 

dW = 1^(4^4 + |^ + ^A3)fc-^/V/^ (41) 

Vr = -leil-A, + A, + 2A,)k-'/'x'/' (42) 

K = -^{28A, + lA, + 5As)k-'/'x'/' (43) 

06 6 
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$ = -x^<^s + — 44 

O X 

^ = -$ - 2eA4A;"^/23,3/2 _ ^^^^ 

We use these results as initial conditions for the system (|T^) to ([T^). Eqs. (|^) and 
(^Tj) show, that the perturbations are in general non-adiabatic. 

Due to the conservation equations (^ and (|^), the constants Ai are not inde- 
pendent. Taking the sum of the ensemble averages of the conservation equation Eq. 
(ID multiplied by /* (the complex conjugate of /p), and its complex conjugate, we 
obtain 



|(i/pp) + ^'(/;/. + /p/:)+2^ 



+3-(/;/p + /p/;) = 0. (46) 

(jj 

Let us discuss the above equation on super-horizon scales where we can neglect the 
second term. Since {d/dt){\fp\'^) = — (|/pP)/t, we find that the real part of {fpfp) is 
negative. Furthermore, Schwarz inequality leads to 

1 1 

lAol > -lAil in the radiation era, and lAol > -lAil in the matter era. (47) 

Similarly, from Eq. (^) we conclude that the real part of {f*fp) is positive and 
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1^2! >- 1^3 1 in the radiation era, and |v42| > -jAsI in the matter era. (48) 



The equality sign is valid, if and only if 



{f;fp + fpp = -2v/(i/pP)(i/pP) • 

We call this condition, which requires that fp and fp are in perfect phase correlation 
"perfect coherence" between fp and fp. On super-horizon scales the spectrum of k'^fy 
is much smaller than the spectrum of fp/t; and thus k'^fv{{fp/t, almost everywhere 
in the space of realizations. Hence, on super-horizon scales energy conservation (Eq. 
I) yields 

fpit) = ^fpiUn) - 4t r ^it')fpit')dt' ■ 

a[t) a[t) Jti„ 

The question of coherence between fp and fp is thus reduced to the question of 
unequal time coherence of fp and fp themselves. Similarly, the coherence between fy 
and fp reduces to the unequal time coherence of each of these functions. We believe 
that for scaling sources, unequal time coherence is reasonably well maintained on 
super-horizon scales and therefore the equal signs in Eqs. (^7]) and ( ^S]) are probably 
valid on sufficiently large scales. Numerical simulations for global scalar fields and 
the large A^ limit (see |TB[) support this hypothesis. 



We now address the effect of unequal time coherence, {fi{t)fj{t')), on the resulting 
power spectrum Ci. To simplify the relevant equations, we neglect here the short 
matter dominated period before decoupling and we also neglect baryons, such that 
c^ = Wr = 1/3. The dynamical components during the tight coupling epoch are thus 
reduced to radiation and seeds. A WKB solution of the evolution equations then 
gives 

D^;\k,x) = -^ r dx'\^x')-^ix')]s^n{{x-x')/^/3) 
V3 "'0 

Vr{k,x) = - r dx'[<^{x')-'^{x)]cos{{x-x)/V3) , (49) 



where we set $(0) = ^(0) = D'^J^ = 0. 

This actually just reformulates our simplified system of equations in terms of two 
integral equations, since the Bardeen potentials $ and \E' are given in terms of 
Dg^\ Vr and the source functions as follows 

$_^ = ^ (^Dir)^^v)A.2e( — f+ — f+f) 

6 + x24^ ^ X ''^^ ^k^i6 + x^)^'^ kiG + x^)^"^^''^ 

'>W + -V;,) + -^$. + 26/_ (50) 



6 + x2'2 ^ X ' 6 + x2 

As we have seen earlier, it is a very bad approximation to replace $ and \1/ by 
the corresponding source potentials $s and "^s, since this does not take care of the 
compensation. 

Applying the Hu and Sugiyama formalism for topological defects [^, we obtain 
within our approximation {Qb = 0, purely radiation dominated) 

-JiiXo - Xdec) C0S( 7W—)f 
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+ [^{Xdec) - ^{Xdec)]je{Xo - Xdec) 
i-xo 

+ / dx[^'{x) - $'(x)] ji{xo - x) 



(51) 



where Xdec = ktdec, Xq = Hq. The times t^ec and to denote the time of decouphng and 
today respectively. The somewhat more involved formula which takes into account 
the presence of baryons and CDM, can be found in [^. 

Assuming coherence, the power spectrum Ci can be calculated by squaring 
A^/(2£ + 1) obtained from Eq. (|5TD, where each variable is replaced by the square 
root of its power spectrum. On the other hand, for totally incoherent perturbations, 
one should instead use IQ 



Cf 



dkk^ 



•^de 



dxPr.{<^-m){ 



Ji[Xo - Xde 



V3 



■ sm 



' Xdec X , 



V3 



-Ji{Xo - Xdec) C0S( 



Xde 



X 



V3 



)y 



+ Jdkk^[P{^ - <!>){tdec)] jKxO - Xdec) 

+ I dkk^[j[{Xo - Xdec)je{Xo - Xdec)[Pr{'^ " ^)]{Xdec)] 

+ fdk r dxk'^Pri^' - $')j£(a;o - x) , 



(52) 



where P{X) = (|Xp)(A;,t) denotes the power spectrum of the variable X and Pr{X) 
is the power spectrum of X integrated over a short time period At (see |^). 

We now want to illustrate the difference of the two approaches in our simplified 
pure radiation model. For pure radiation we can derive the following second order 
equation for -D^''^: 



3 (6 + x^)x ' 6 + x2 



with general solution 

D^;'^ = -e r f(x')G(x,x')dx' 
^ 3 Jo 



ni;) 



X 



k^ie + x^)^'^ 



3x 



3 

8. „ , x' 



k(6 + X 



2\Jv~t~Jtt 



$. 



(53) 



(54) 



where /(x) = firix) + (x^/6e)$s/ and G{x,x') denotes the Green's function 



G{x,x') 



VSx' 



1^12 + xx') sin 



(6 + x"^)x 
The power spectrum of D)^' is therefore 
64 



x — x 



V3 



2v3(x + x) cos 






v^ 



.(55) 



a 
ex rx 



(|dM|2)= 2 / / dx'dx"{f{x')nx"))G{x,x')G{x,x") 
^ 9 Jo Jo 

Assuming total coherence, Eq. (^) takes the form 

(I^^H = [lejyx'^{\f{x')\^)G{x,x')]\ 



(56) 



(57) 
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On the other hand, assuming complete decoherence, 

rx+h.x 



(/(x)r(x')) = 6{x - x') / dx{\f{x)\') , (58) 

Jx 

leads to the power spectrum 

(I^WH = -ehj^ dx\\f{x')\')G\xy) , (59) 

where we have chosen Ax = rx. We further assume also complete decoherence 
between different source functions, 

(/./.) = (/./p) = {Up) = . (60) 

In Figs, la and lb we plot \D^g'^\'^k^ versus ktdec under the assumption of total coher- 
ence and complete decoherence respectively. The role of the coherence assumption 
on the characteristics of the power spectrum is shown in Fig. 2. Clearly, complete 
decoherence shifts the first acoustic peak to smaller angular scales, and reduces 
substantially its height. Furthermore, secondary peaks are completely washed out. 
A realistic defect model will always lay somewhere between these two extremes. 
We suppose however, motivated by numerical simulations of textures and the large 
A^ limit, that the texture example is closer to the completely coherent case. In the 
next section we thus restrict ourselves to perfect coherence. 



4 Numerical examples 



In this section we study how the characteristics of the acoustic peaks depend on the 
values of the dimensionless constants Ai and the form of the functions Fj, which 
determine the power spectra of the seed functions (see Eq. (PUf)). A crucial question 
is whether there is a set of parameters for which the position and amplitude of the 
primary acoustic peak are similar to those predicted by an adiabatic inflationary 
model. 

As we discussed earlier, the functions Fi are normalized such that -Fj(O) = 1, and 

Fi —>■ ioT X ^ oo. Numerical simulations for global textures [jlOl suggest that in 
the case of global scalar fields, the functions Fi have power law decay. However, for 
generic scaling sources, one could also consider the case of exponential decay. As we 
shall show, the form of these functions affects the features of the power spectrum 
significantly. In general, we find that if Fi have an exponential decay, the position 
of the primary peak is within the range predicted by adiabatic inflationary models, 
at £ ~ 220. On the other hand, if Fi have a power law decay, as it seems to be for 
global topological defects, the position of the first acoustic peak is clearly shifted 
to smaller angular scales, at around i ~ 300 to 400. The amplitude of fluctuations 
decreases by up to a factor of 500 if we choose an exponential decay law for the seed 
functions. This is due to the fact that in this case the decay is very fast and erases 
almost all substantial seed contribution. The sensitivity of the overall amplitude 
of CMB perturbations on the parameters is extremely important especially if one 
wants to rule out defect models with biasing arguments! 

In figures 3 to 8, we show the resulting power spectra for different set of dimen- 
sionless constants Ai, with the seed functions having either exponential or power law 
decay. The dashed line indicates the SW contribution, the dashed-dotted line is the 
contribution from the acoustic peaks, and their sum is drawn as solid line. 
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In Fig. 3a we choose an exponential decay for all the seed functions, Fi{x) = 
exp(— x^), and a set of constants Ai = 3, ^3 = —0.6,^4 = 0. The position of the 
primary peak is at £ ~ 200, while the relevant height of the first peak with respect 
to the SW plateau is ~ 25. Using the same set of constants A^, however choosing 
a power law decay for the seed functions, -Fj(x) = [1. + (l./(27r)^)x^]~"'^, we see in 
Fig. 3b that the peak is displaced to smaller angular scales, at ~ 330, while the 
relative amplitude of the acoustic peak with respect to the SW plateau remains 
~ 25. Also the features of the secondary peaks are different. While in Fig. 3a the 
second and third peaks have almost the same height, in Fig. 3b the second peak has 
almost completely disappeared. In both cases, the spectral index of the plateau, in 
the range £~2 — 20isn~l, consistent with observations. 

Now, selecting a slightly different set of dimensionless constants, we see that the 
predicted power spectrum is very different. In Fig. 4 we show the power spectrum 
for Ai = 3, ^3 = —0.7, ^4 = 0, and the same power law decay for the functions Fj 
as in Fig. 3b. The primary peak is again at rather large angular scales, ~ 350, but 
the height of the peak is different, and the spectral index clearly deviates from 1. 

A very interesting case is displayed in Fig. 5a, where we see that both the position 
and the amplitude of the first acoustic peak, agreed with those predicted by a generic 
inflationary model. Here, Ai = 3, ^3 = 1,^4 = 2, and Fi{x) = exp(— x^). The 
primary peak is at £ ~ 200 and its relative amplitude is at ~ 4. The second 
peak is almost completely washed out and the spectral index in the range i ~ 
2 — 20 is very close to 1. This power spectrum, where perturbations are generated 
by scaling seeds, is quite similar to one resulting from an adiabatic inflationary 
model. Considering the same set of parameters and a somewhat slower exponential 
decay for the functions F, given by Fj(x) = exp{—0.5x'^), we find (Fig. 5b) that 
both, the position and relative amplitude of the first peak with respect to the SW 
plateau, remain the same as in Fig. 5a, whereas the Sachs Wolfe plateau is somewhat 
prolonged. This simple example shows that it may well be possible to "manufacture" 
inflationary spectra by a suitable choice of seed functions. A point which has already 
been realized in Ref. |TP[. It is thus extremely important to further constrain the 



seed functions of defect models by numerical simulations and/or the large A^ limit. It 
may well be that the requirement of power law decay of the seed functions, excludes 
the inflationary position of the first acoustic peak. With the same parameters Ai 
and a power law decay Fi{x) = [1 + (l/(27r)^)s^]~-'^, the primary peak is at £ ~ 320, 
while the relative height of the peak is about 8, Fig. 5c. 

We have also considered the values for the parameters Ai which are suggested 
by the conservation equations and perfect coherence, choosing a power law decay, 
Fi{x) = [1 + (l/(27r)^)a;^]"^ We find that the sign of ^3 does not affect the features 
of the power spectrum. In Fig. 6a Ai = 3, ^3 = —1/6,A4 = l/(2(27r)^), while in 
Fig. 6b Ai = 3, A3 = 1/6,^4 = l/(2(27r)^). In both these cases, the peak is at 
i ~ 300, and its amplitude is ~ 8. 

Finally, to illustrate the variety of results which can be obtained by parameter 
variation within a simple family of seed functions, we show a rather extreme case 
in Fig. 7, where Ai = 3,^3 = — l.,A4 = and Fi{x) = exp{—x'^). Here we see no 
acoustic peaks at all. However, this is a rather particular case, since with this choice 
of parameters, all the variables $s, \l/s, D^J'\ D^^\ V^ and Vc vanish initially (see Eqs. 



(HI) to (ID). 

To analyze the dependence of the characteristics of the power spectrum on the 
seed functions in a somewhat more systematic way, we have calculated the C^'s 
for a grid of values — Ai < ^3, A^ < Ai with spacing 0.2y4i and fixed functions 
Fi = 1/(1 + (x/(27r))2). We fitted the resulting Q for 2 < £ < 20 to the simple power 
law behaviour arising in inflationary models, Ce oc T{i+ {n — l)/2)/r(£+ (5 — n)/2). 
We find that for IA4I > 0.4y4i, the x^ of the fit is unacceptably high: x^ ~ 3 to 4 for 
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|y44| = 0.4yli, and more than 10 for even larger anisotropic stresses. We allow for 
a relative error of 0.05. However, x^ depends only weakly on the value of A^ (see 
Fig. 8). We therefore restrict the parameter range for A4 to —O.SAi < A4 < O.SAi. 
The spectral index is in good agreement with observations, I < n < 1.4 (see Figs. 9a 
and 9b). 

We find positions of the first acoustic peak in the range 260 < ipeak < 520. For 
the choice of seed functions with power law decay on sub-horizon scales (which is also 
indicated from numerical simulations and from the large N limit), we never obtain 
the peak at the adiabatic inflationary position of ~ 220, and values Ipeak < 300 are 
only found for very small A4 (see Fig. 10). 

We define the quantity ipeaki^peak + 1)C*£ eak/i^^^^w) as a measure for the height 
of the acoustic peak. This quantity is very model dependent and assumes, within 
the small class of models investigated in our parameter study, all values between 0.1 
(for ^4 = and ^43 = ^41, i.e., virtually no discernible peak) and 11 (see Figs. 11a 
and lib). For fixed A3, the peak height is a steeply raising function of {A^l. Only 
values IA4I < 0.1 Ai lead to peak heights below 6. 

It is interesting to note, that also the absolute amplitude of the spectrum is 
sensitively depending on the ratios ^4/^41 and ^43/^1. For ^4 = the amplitude 
llOCio varies from 0m2e^Aj for A3 = -0.2Ai to O.SeM^ for A3 = A^. If A^ ^ 0, 
the amplitude does not depend very strongly on ^3 and grows from ~ OAe^Af for 
\A^\ = O.lAi to ~ Ae^Aj for \A^\ = 0.3Ai (see fig. 12). This finding is important 
for the biasing problem of structure formation. Sometimes, defect models have been 
claimed to be ruled out, since they would not lead to large enough matter density 
fluctuations, if normalized to the COBE experiment on very large scales. This 
normalization fixes the only free parameter of a given model, namely the symmetry 
breaking scale and therefore e. In our work, we have seen that the Sachs Wolfe 
fluctuations in the CMB are largely governed by A4, the amplitude of anisotropic 

stresses. The density fluctuations in the dark matter, however, are induced by 0^ 
alone (see, e.g. [|r^), which is determined entirely by Ai and A2. Defect models, 
with somewhat small anisotropic stresses, e.g., A4 < 0.05Ai, which are actually quite 
natural, but difficult to resolve numerically, may explain the different bias factors 
obtained from numerical simulations in lEOl, p!0|. 



5 Conclusions 



In this paper we analyzed with some generality the CMB anisotropics induced in 
models with scaling sources for an f2 = 1, cold dark matter cosmological model. 

Within the framework of gauge invariant perturbation theory it turns out that 
compensation is a consequence of "natural" initial conditions. By "natural" we mean 
that we only consider that part of the solution induced by the source itself and do 
not add an arbitrary homogeneous contribution. In this case, we have found that 
the total Bardeen potentials are reduced by a factor x^ = {ktY with respect to the 
potentials generated by the source alone. One may think at first sight that such a 
result is unphysical, acausal, however it just refiects that also the initial condition 
of a perfect Friedmann universe is acausal. 

Even restricting ourselves to the case of scaling sources, we found that the result- 
ing power spectrum depends significantly on the model parameters. In particular, 
if the seed functions decay exponentially, the position of the first acoustic peak is at 
£ ~ 220 as in infiationary models. Adjusting the amplitude of the seed functions. 
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we can also obtain a peak height consistent with inflationary perturbations. 

On the other hand, if the seed functions have a power law decay, as numerical 



simulations of global textures |]T0| as well as the large A^ limit [|T8] indicate, the 
position of the first peak is witmn the range 260 < ipeak < 500. Its amplitude 
depends sensitively on the parameters of the seed functions which, for a specific 
model, have to be determined by involved numerical simulations. In our analysis we 

encountered amplitudes in the range 0.1 < ipeaki^peak + 1)C|'^" /(llOCio) < 25. 

We also found that the total amplitude of CMB anisotropics produced depends 
strongly on the amplitude of anisotropic stresses of the seed. Whereas, the source 
term leading to CDM density fluctuations is given by fp + 3/p, i.e, determined by 
Ai + 3A2 on large scales. Therefore, the relation between the COBE normalization 
of the model and the bias factor depends sensitively on the ratio A^/i^Ai + 3A2), 
which may depend on details of the model. 

These results are obtained under the assumption of perfect coherence. This 
hypothesis seems reasonable for global scalar fields as also the large A^ limit indicates 
|lq| . On the other hand, assuming complete decoherence, the position of the first 
peak is shifted to smaller angular scales, its amplitude is reduced and secondary 
peaks are washed out. A realistic situation may lay somewhere between the two 
extremes. 

The examples with power law seed functions which we discussed in this paper 
were motivated by numerical simulations of global textures with vacuum manifold 
S^. Apart from the scaling behaviour on very large and very small scales, which 
should be the same for all global defects, we do not know to what extend the seed 
functions depend on this particular choice. 
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APPENDIX 

A Compensation for relativistic collisionless par- 
ticles 



In this appendix we show how compensation arises in the case of relativistic colli- 
sionless particles. 

We consider a universe dominated by massless [i.e. relativistic) collisionless 
particles, with scalar perturbations induced by seeds. We assume seeds consisting of 
massless particles conformally coupled to gravity. In this case, the time dependence 
of the seed functions is given by /, oc 1/a^ oc 1/t^ {t denotes conformal time). 

The evolution of perturbations is determined by the collisionless Bolztmann equa- 
tion, which reads llT2i 



dtM{fx,k,t)+ikfxM = ikfx[^ - '^]{k) , (Al) 
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where k denotes the wave number and /i = n ■ k/fc; n stands for the momentum 
direction of the relativistic particles. A^ is a gauge-invariant perturbation variable 
for the energy integrated one-particle distribution function, 



M 



TT 



P 



p^dp6f 



T'' = J^f{p)p'f 



Using the general definition of the energy momentum tensor, 

we obtain 

Dg = 2 f Mdfi 

3i /■! 

V = — Mfxdfx 

2 j-i 

n = 3 /" (1 - 3fi^)Mdfi . 
The gravitational perturbation equations on super-horizon scales yield 
$ 



4x6 

-$ - 2eU - 4n 



x^ 



(A2) 
(A3) 
(A4) 

(A5) 
(A6) 



where ^s = ^k ^[/p + (3/t)/t,]. Furthermore, the conservation equation, D' 



— (4/3)1^, tells us that Dg oc xV , so that we may neglect the Dg term in Eq. ([A5| 



Inserting Eqs. ([A5|, ^3) in Eq. (|A1D, we obtain 



dyM{y,x) + iM =i[ — + -x^^, + 2e/^ + — ](fc) , 

X 3 x^ 

where y = jix. 

For X -C 1, and thus y ^ 1, we make the ansatz (see 



^ 



M = x^[ciy + C2y'' + 0{y')] , 

where ci, C2 and (3 are constants. Inserting this ansatz in Eqs. 
we get 



(AT) 



(A8) 
D, (^) and (^, 



D 



9 

n 



-C2X 



/3+2 



2CiX' 



/3+1 



-(/3 + 2)c2X 



/3+1 



--CoX 

5 



/3+2 



(A9) 
(AlO) 
(All) 



where the second equality in Eq. ( |A10|) is obtained from the energy conservation 
equation (the zeroth moment of Eq. ( |A1|) ) . The first moment of Eq. ( [Al|) (momen- 
tum conservation) implies 



[(/5 + 2)(/? + 3) + -]C2X^ = -x2$, + 2e^ 



(A12) 
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lieU ~ a;2$„ Eq. (^) leads to 0{V/x) = 0{U/x^) = 0{x^^s) = 0{x^^s). Thus, 
(9($) = (9(x^$s) and (9(\1/) = 0{x'^'^s) which means that we find compensation. In 
the other case, /^ ^ x^$s, the right hand side of Eq. ( |A12| ) behaves hke e/^ oc 1/x^, 
implying [3 = —2. Eq. ( |A12| ) then leads to 



-li/x^ = \c2X^ ~ 2e/^ and V ~ . 





Inserting this in Eqs. ( |A5| , [A 61 ) yields 

^ = -$ = -ia;2<l>. , (A13) 



and thus (9($) = (9(a;^$s) and 0{^) < 0{x'^'^s), leading again to compensation. 

We thus have shown that in this example, compensation is present even if the 
anisotropic stresses of the seeds are not suppressed. 
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Figure 1: The acoustic fluctuations in the photon density spectrum are shown for 
the case of perfect coherence, (top) and complete decoherence (bottom). 



19 



100 




Figure 2: The resulting spectrum of CMB anisotropies from the photon density per- 
turbations given in Fig. 1 for perfect coherence (sohd hne) and complete decoherence 
(dashed line). 
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Figure 3: Here and in the subsequent figures, the CMB anisotropies are shown in 
units of e^ = {AttGi]'^)'^. The Sachs Wolfe contribution alone is indicated by a dashed 
line and the coherent sum of acoustic and Doppler terms are shown as dot-dashed 
curve. The solid line is the incoherent sum of these two contributions. Here the seed 
functions are determined by the choice Ai = 3, A^ = —0.6 and A^ = 0. In the top 
frame the seed functions decay exponentially, while they decay like a power law in 
the bottom one. 
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Figure 4: As Fig. 3b, but with A3 = —0.7. In this regime (A4 = 0, A3 ~ —(1/4 to 
l/3)Ai) the resulting spectrum depends very sensitively on A3. While A3 = —0.6 
leads to a perfectly reasonable spectrum with a somewhat high first acoustic peak, 
this spectrum is excluded by observations due to its "bump" in the Sachs Wolfe 
plateau and the absence of a distinctive acoustic peak. 
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Figure 5: The CMB anisotropy spectrum for parameters Ai = 3, A3 = 1, A4 = 2 
and exponential decay, Fi oc exp(— a;^) (top) and Fi oc exp(— a;^/2) (middle). The 
corresponding spectrum for seed functions with power law decay is shown in the 
bottom frame. The position and relative amplitude of the first acoustic peak of 
the spectra (top) and (middle) is comp^ble with an inflationary spectrum. This 
simple exaniDle hints that it mav be possible to "manufacture" inflationary spectra 
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Figure 6: The anisotropy spectra for Ai = 3, A4 = Ai/6(27r)^ and A^ = Ai/18 
(top) respectively A3 = —Ai/18 (bottom) are shown. The seed functions are chosen 
to have power law decay. These are the values Ai which can be inferred from energy 
momentum conservation under the assumption of perfect coherence (see text). The 
sign of A3 cannot be deduced, but we see that the results do not depend on it. 
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Figure 7: Even extremely strange spectra, like this one with a negative spectral 
index and without acoustic peaks can be obtained. For this result we have chosen 
exponentially decaying source functions, Fi ex exp(— x^) and Ai = 3, A3 = —1, 

^4 = 0. 
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Figure 8: The CMB anisotropies on large scales, i < 20 from 122 models with 
scaling sources with structure function amplitudes (for details see text) in the regime 
—1 < A4/A1 < 1 and —1 < A^/Ai < 1 are fitted to simple power law spectra with 
spectral index —0.5 < n < 2.5. The x^ of the fit (allowing for 5% relative error) is 
shown as function of A4. 
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Figure 9: The spectral index for different models with —0.3 < A4/A1 < 0.3 and 
— 1 < A3/A1 < 1 is shown as function of A4 (top). The fact that there are less than 
six different circles visible for some values of A4 is due to the discrete spacing of 
about 0.06 in n. The bottom frame shows the same results as a function of A^. 
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Figure 10: The position of the first acoustic peak is shown as function of A4/A1 for 
different values of A3 in the range —Ai < A3 < Ai. 
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Figure 11: The amplitude of the first acoustic peak over the Sachs Wolfe plateau 
is shown as function of A4/A1 (top) for different values of A3 in the range —Ai < 
A3 < Ai. The bottom frame shows the same results as a function of A3. 
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Figure 12: The amplitude llOCio is given (in units of e^Af) as a function of the 
amphtude of anisotropic stresses, A4/A1. It varies over about 3 orders of magnitudes 
and can become substantially smaller than 1, especially for very small values of A4. 
The significance of this finding for the biasing problem is discussed in the text. 
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